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scheme.  This  work  can  form  a  basis  to  unravel  the  underlying  science  of  fly-ear  mechanisms  and  provide  a  framework  for  the  development 
of  bio-inspired  localization,  tracking,  and  navigation  systems. 
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I.  Introduction 

Animal  hearing  organs  are  governed  by  a  fundamental  size  constraint;  the  smaller  the 
organ  size,  the  smaller  the  available  directional  cues  for  directional  hearing.  However,  with  an 
auditory  organ  separation  of  only  520  pm,  it  is  remarkable  that  the  parasitic  fly  Ormia  ochracea 
can  achieve  a  human-like  localization  precision.  The  key  to  this  fly’s  phenomenal  directional 
hearing  ability  is  believed  to  be  due  to  the  mechanical  coupling  between  its  two  eardrums.  The 
overall  goal  of  this  AFOSR  sponsored  research  program  is  to  achieve  a  fundamental  understanding  of  the 
mechanism  the  fly  ear  and  use  this  understanding  to  develop  fly-ear  inspired  miniature  acoustic  sensors 
for  target  detection  and  localization.  Over  the  three-year  research  program,  significant  progresses  have 
been  made  toward  achieving  this  goal.  First,  through  mechanics  and  optimization  analysis,  we  have 
discovered  an  important  biological  conclusion:  the  fly  ear  can  be  viewed  as  a  nature-designed  optimal 
structure  that  is  endowed  with  the  dual  optimality  characteristic:  maximum  average  directional  sensitivity 
and  minimum  nonlinearity  over  the  azimuth  range  of  -30  degrees  to  30  degrees  achieved  simultaneously 
at  the  fly  ear’s  working  frequency  of  5  kHz.  Furthermore,  a  micro-scale  fly-ear  inspired  directional 
microphone  has  been  developed  to  fully  capture  the  dual  optimality  characteristic  of  the  fly  ear.  A  lab-on- 
a-chip  optical  detection  system  has  also  been  developed  to  realize  a  fully  integrated  and  miniaturized 
sensor  system.  Moreover,  a  localization  strategy  inspired  by  the  fly’s  localization/lateralization  scheme 
has  been  developed  and  employed  experimentally.  In  addition,  it  has  also  been  found  that  not  only  the  fly 
ear  represents  a  nature  designed  optimal  structure;  the  fly’s  localization  scheme  may  also  resemble  an 
optimal  control  scheme.  This  work  can  form  a  basis  to  unravel  the  underlying  science  of  fly-ear 
mechanisms  and  provide  a  framework  for  the  development  of  bio-inspired  localization,  tracking,  and 
navigation  systems. 

II.  Technical  Summary  of  Accomplishments 

A.  Fundamental  understanding  of  the  fly’s  dual  optimality  characteristics 

In  this  project,  we  have  furthered  the  understanding  of  fly-ear  mechanisms  and  used  this 
understanding  to  develop  fly-ear  inspired  sensors.  We  have  found  the  answers  to  the  following 
fundamental  questions:  i)  how  have  the  structural  parameters  of  the  fly  ears  been  tailored  for  achieving  its 
superior  localization  ability  at  5  kHz,  ii)  does  the  fly  ear  represent  a  nature-designed  optimal  structure  to 
facilitate  the  fly’s  unique  localization-lateralization  scheme,  and  iii)  can  a  synthetic  device  be  developed 
to  truly  replicate  the  fly-ear  characteristics. 

The  mechanics  model  used  in  this  study  is  an  equivalent  two-degree-of-freedom  model  along 
with  its  relevant  parameters  (Miles  et  al,  1 995 ’).  As  shown  in  Figure  1,  the  two  tympana  are  modeled  as 
two  rigid  bars  connected  by  a  torsion  spring  k 3.  The  outer  end  of  each  bar  is  supported  by  a  translational 
spring  k]  or  k2  ,  which  is  equivalent  to  the  tympanum  stiffness.  Dashpots  c\,  c2  and  c3  are  added  to  account 
for  the  damping  ratios  of  the  system.  All  the  parameters  used  to  study  the  fly  ear  structure  are  listed  in 
Table  1. 
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Figure  1:  Equivalent  2-DOF  model  of  the  fly  ear. 


The  governing  equations  of  the  coupled  system  can  be  obtained  as  follows: 
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Table  1:  Parameters  used  in  the  2-DOF  model  for  the  fly  ear  (Miles,  et  al.,  1995) 


Parameters 

Values 

Mass  of  bar  m 

2.88xlO'10  kg 

Translational  spring  kuk.2 

0.576  N/m 

Translation  dashpot  c/,c2 

1 . 1 5 x  10"5  Ns/m 

Torsional  spring  ks 

5.18  N/m 

Torsional  dashpot  cq 

2.88xl0'5  Ns/m 

Separation  of  force  locations  d 

1.2x1  O'3  m 

Tympanum  area  s 

0.288xl0'6  m2 

Excitation  frequency  ax 

3.14xl04  rad/s  (5  kHz) 

Sound  speed  c 

344  m/s 

In  the  case  of  free  vibration,  the  natural  frequencies  and  mode  shapes  of  the  system  can  be  found 


as: 


OX 


=  yjkjm,  co2  =  ^{ky  +2k^)/m 

Vi=[-i  l]r,v2=[l  if, 


where  the  two  diaphragms  move  out  of  phase  in  the  first  mode  (rocking  mode),  and  in  phase 
mode  (bending  mode). 


(3) 
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in  the  second 


Following  an  approach  different  from  that  used  in  the  existing  literature,  modal  analysis  is  carried 
out  and  this  allows  us  to  obtain  analytical  solutions  for  the  directional  cues  of  the  fly  ear  and  understand 
how  the  structural  parameters  of  the  fly  ear  affect  its  performance.  By  using  modal  analysis2,  the 
responses  to  an  acoustic  stimulus  p0ejai>  (po  is  the  pressure  reference  on  the  midline)  are  obtained  as 


Jxi  1  =  \Al \e  jat  =  \T  +  J  tan(^/2)| 

jx2j  jH2  j  1-Q2  +/2|jQ  jr  -  _/  tan(^/2)| 
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Here,  T  represents  the  relative  contribution  of  the  two  modes  subject  to  unit  modal  force,  which  is 
dependent  on  the  natural  frequency  ratio  rj,  the  damping  factors  §ij2,  and  the  normalized  excitation 


frequency  £2.  (p  represents  the  phase  difference  of  the  incident  sound  field  applied  to  the  two  mass-spring 
systems. 

Then,  mechanical  interaural  intensity  difference  (mllD)  and  mechanical  interaural  phase 
difference  (mIPD)  between  the  two  eardrums  can  be  analytically  obtained  as 


mllD  =  20  log10 


(8) 
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Equivalent  to  mlTD,  the  mechanical  interaural  phase  difference  (mIPD)  is  chosen  as  the 
directional  cue  for  investigation,  since  it  is  a  dimensionless  measure  that  is  independent  of  the  sound 
wavelength  and  sound  speed. 

On  the  midline  (0  =  0),  both  mllD  and  mIPD  are  equal  to  zero.  In  our  study,  mIPD  is  chosen  as 
the  main  directional  cue,  for  which  the  directional  sensitivity  (DS)  is  defined  as  the  first  derivative  of 
mIPD  with  respect  to  azimuth: 


DS  =  dmIPD/dd  ■ 


(10) 


It  is  worth  noting  that  the  localization  performance  does  not  only  depend  on  the  absolute  value  of 
mIPD,  but  more  importantly,  it  is  the  directional  sensitivity,  that  determines  how  accurately  the  fly  can 
pinpoint  a  source. 

Based  on  the  analytical  solutions  of  mIPD  and  DS,  it  is  easy  to  see  that  mIPD  and  DS  solely  depend 
on  two  parameters:  the  modal  response  ratio  r  and  the  initial  phase  difference  (j)  .  The  former  relates  to 
the  resonance  frequency  ratio  17  (the  ratio  of  the  bending  mode  natural  frequency  <w2  to  the  rocking  mode 
natural  frequency  o.)\),  normalized  working  frequency  Q  (the  ratio  of  the  sound-source  frequency  to  the 
rocking  mode  natural  frequency),  and  damping  ratios  £1  and  §2.  The  latter  relates  to  the  separation-to- 
wavelength  ratio  %  (the  ratio  of  the  membrane  center  separation  d  to  the  sound-source  wavelength  A)  and 
the  sound-source  azimuth  6. 

By  looking  at  the  directional  sensitivity,  as  shown  in  Figure  2,  we  can  find  out  that  in  the  vicinity 
of  midline  (-3O°<0<3O°),  DS  is  not  only  very  flat  at  the  fly’s  working  frequency  5kHz,  but  it  is  higher 
than  what  is  achievable  at  any  other  frequency.  In  other  words,  if  we  use  a  straight  line  to  approximate 
mIPD  in  the  azimuth  range  -3O°<0<3O°,  the  linearity  is  best  at  5kHz,  and  the  slope  is  maximal.  The 
localization/lateralization  scheme  takes  full  advantage  of  the  linear  range  so  that  the  fly  can  accurately 
pinpoint  the  crickets. 


Figure  2:  Directional  sensitivity  (DS)  as  a  function  of  azimuth  for  various  frequencies.  5kHz  is  the 
calling  song  frequency  of  the  fly’s  host  crickets. 


Further,  two  performance  metrics  for  further  investigation  can  be  defined:  the  nonlinearity  (NL) 
of  mIPD  and  the  average  DS  (ADS)  over  the  azimuth  range  of  -30°  to  30°.  When  these  two  metrics  of  the 
fly  ear  are  plotted  in  the  frequency  domain  in  Figure  3,  we  are  able  to  observe  an  interesting  result:  the 
minimum  NL  and  the  maximum  ADS  are  achieved  simultaneously  at  the  frequency  of  5kHz.  This  result 
suggests  that  the  fly  ear  is  endowed  with  a  dual  optimality  characteristic  at  its  working  frequency. 
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Figure  3:  Directional  sensitivity  and  nonlinearity  of  the  fly  ear  in  the  frequency  domain. 


We  further  explore  as  to  how  the  structural  parameters  of  the  fly  ear  are  tailored  to  achieve  such  a 
dual  optimality  characteristic  and  whether  a  synthetic  device  endowed  with  the  fly  ear’s  dual  optimality 
characteristic  can  be  developed.  An  optimal  design  problem  is  formulated  to  seek  solutions  in  the  design 
space,  which  will  meet  the  objective  of  achieving  dual  optimality  (i.e.,  minimal  NL  and  maximal  ADS) 
over  the  azimuth  of  -30°  to  30°  at  the  same  sound  frequency  (i.e.,  the  optimal  working  frequency).  For 
simplicity,  the  problem  is  constrained  in  one  dimension  (azimuth).  Note  that  there  are  several  key 
parameters,  the  resonance  frequency  ratio  r],  the  separation-to-wavelength  ratio  x,  and  the  damping  ratios 
that  can  influence  the  NL  and  ADS.  In  Figure  4,  the  rocking  and  bending  mode  natural  frequencies  that 
ensures  the  dual  optimality  characteristic  are  plotted  as  a  function  of  separation-to-wavelength  ratio  %  for 
two  damping  scenarios:  §i=0.89  and  §2=1.23  for  the  fly  ear  case,  and  §i  =0.17  and  §2  =0.05  for  the  other 
case.  Given  the  fly-ear’s  geometric  dimensions  and  its  working  frequency  of  5kHz,  following  the  two 
curves  in  Figure  4,  the  natural  frequencies  (rocking  mode:  6.99kHz,  bending  mode:  30.1kHz)  required  to 
achieve  the  dual  optimality  can  be  obtained,  which  are  found  to  be  consistent  with  the  experimental  data 
reported  in  the  literature  (7.12kHz  and  31.0kHz).  This  finding  provides  the  evidence  for  the  statement  that 
the  fly  ear  indeed  represents  a  nature-designed  optimal  structure  that  can  simultaneously  achieve  the 
maximum  DS  and  the  minimum  NL  at  its  working  frequency  of  5kHz. 

Furthermore,  to  achieve  such  a  dual  optimality  characteristic,  proper  contributions  from  both  the 
rocking  and  the  bending  modes  (i.e.,  resonance  ratio  77=4 . 3 6)  are  necessary.  It  is  noted  that  the  resonance 
ratio  rj  is  related  the  stiffness  ratio  cr  =  kj/k  1  by  J72=l+2cr,  which  is  the  key  non-dimensional  structural 
parameter  that  determines  how  strongly  the  two  ear  membranes  are  coupled  (i.e.,  the  coupling  strength). 
As  shown  in  of  Figure  5,  given  the  same  separation  to  wavelength  ratio  as  the  fly  ear,  if  the  coupling  is 
weak  (e.g.,  r]=2),  the  amplification  of  phase  difference  is  not  significant,  i.e.,  DS  is  too  small  to  achieve 
the  maximal  value.  On  the  other  hand,  when  the  coupling  is  rigid  (e.g,  rj=20),  even  though  the  mIPD  can 
be  significantly  amplified,  it  saturates  at  180°  when  6  is  slightly  off  the  midline,  and  thus,  it  is  not 
possible  to  distinguish  the  azimuths  and  the  minimum  NL  is  not  achievable.  Only  when  the  coupling  is 
“medium”  36),  the  fly  ear  can  achieve  a  balance  between  DS  and  NL,  and  thus,  the  dual  optimality 
can  be  achieved  at  its  working  frequency.  This  suggests  that  the  structural  parameters  of  the  fly  ear  may 


have  been  well  adapted  in  the  course  of  evolution,  resulting  in  an  appropriate  coupling  strength  for 
achieving  the  dual  optimality  characteristic. 


Wavelength-to-separation  ratio  X/d 

Figure  4:  Design  curves  to  achieve  fly  ear’s  dual  optimality.  Two  damping  scenarios  are  shown: 

=  0.89  and  lj2  =  1-23  for  the  high  damping  case,  and  =  0.17  and  f2  =  0.05  for  the  low  damping  case. 
For  any  given  wavelength-to-separation  ratio  in  each  case,  the  two  natural  frequencies  for  the 
rocking  and  bending  modes  can  be  obtained  using  the  optimal  design  curves  here,  e.g.  the  fly  ear  in 
the  high  damping  case,  and  the  synthetic  MEMS  device  in  the  low  damping  case. 
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Figure  5:  Phase  difference  with  respect  to  incident  azimuth  for  different  coupling  strength.  In  order 
to  achieve  dual-optimality,  the  coupling  strength  of  the  bridge  (defined  by  the  ratio  of  two  natural 
frequencies)  needs  to  be  medium,  neither  too  soft  nor  too  rigid. 


The  dual  optimality  characteristic  can  provide  a  basis  for  understanding  the  fly’s  superior 
directional  hearing  capability  as  well  as  its  unique  localization-lateralization  scheme.  First,  although  the 
absolute  value  of  the  mIPD  is  maximal  at  the  two  extreme  positions  (0=±9O°),  the  corresponding  DS  is 
close  to  zero  at  these  positions  and  the  maximal  DS  is  actually  achieved  at  the  vicinity  of  the  midline  (0 
=0°).  Therefore,  the  fly  naturally  turns  the  head  front  (midline  of  the  ear)  towards  the  source  so  that  the 
maximum  DS  (i.e.,  maximal  slope)  can  be  achieved  to  ensure  the  best  localization  precision.  This  is 
similar  to  a  related  finding  reported  for  an  Egypt  fruit  bat,  which  uses  not  the  maximal  sonar  beam 
intensity  but  its  maximal  slope  for  target  localization3.  Second,  mIPD  is  a  linear  function  of  azimuth  in 
the  range  from  -30°  to  30°,  which  is  coincident  with  the  sigmoid  relationship  of  fly’s  turning  speed  with 


respect  to  the  azimuth,  obtained  in  the  phonotactic  experiments  on  the  fly[4].  Given  the  limited  neural 
processing  ability  of  the  fly,  a  linear  and  maximal  DS  can  certainly  help  the  fly  perform  the  localization 
task  more  accurately  and  more  efficiently  for  the  azimuths  from  -30°  to  30°.  Therefore,  in  this  sense,  it  is 
not  only  the  mechanical  coupling  mechanism  that  helps  the  fly  ear  obtain  significantly  amplified 
directional  cues,  but  more  importantly,  the  structural  parameters  of  the  fly  ear  have  been  tailored  to 
achieve  the  dual  optimality  characteristic  at  5kHz,  facilitating  a  unique  localization-lateralization  scheme 
for  accurately  pinpointing  its  host. 

B.  Development  of  Bio-Inspired  Acoustic  Sensor  with  Dual  Optimality  Characteristics 


Inspired  by  the  fly  ear’s  dual  optimality  characteristics,  fly-ear  inspired  sensors  have  been 
developed  and  experimentally  demonstrated.  The  fly-ear  inspired  directional  microphone  consists  of  two 
clamped  circular  membranes  (material:  silicon,  radius  of  500  pm  and  thickness  of  0.5  pm)  and  a  coupling 
bridge  (material:  alternating  SiCF  and  SI3N4  layers,  size:  1250pm  x300pmx3.2pm)  that  is  pivoted  in  the 
middle  and  connects  the  two  membrane  centers.  It  is  noted  that  the  separation  between  diaphragm  centers 
is  designed  to  be  1.25  mm,  same  as  that  of  the  fly  ear.  The  membrane  structure  is  integrated  with  another 
layer  of  silicon,  serving  as  holders  for  two  optical  fibers,  which  are  used  to  detect  the  acoustic  pressure 
induced  membrane  deflection.  The  bio-inspired  MEMS  directional  microphone  uses  a  low  coherence 
fiber  optic  interferometer  (LCFOI)  system  to  detect  the  displacement  of  the  two  membrane  centers.  As 
shown  in  Figure  6,  a  superluminescent  light  emitting  diode  (SLED)  is  used  as  the  broadband  light  source. 
The  sensing  interferometer  is  a  Fabry-Perot  (FP)  cavity  formed  between  each  membrane  and  the 
corresponding  fiber  tip  as  the  sensing  interferometer;  the  read-out  interferometer  is  a  fiber  FP  tunable 
filter  (Micro  Optics  FFP-TF2).  In  order  to  achieve  maximum  sensitivity,  biases  are  applied  to  the  tunable 
filters  to  adjust  the  initial  working  position  to  be  at  quadrature  points.  Compared  with  conventional  fiber 
optic  interferometers,  this  detection  system  can  substantially  reduce  the  phase  noise  associated  with  light 
wavelength  stability.  In  addition,  using  a  much  shorter  cavity  length,  it  can  achieve  high  spatial  resolution 
and  great  reduction  in  temperature  cross  effects. 
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Figure  6:  Low  coherence  fiber  optic  interferometer  system  for  the  mico-fabricated  directional 
microphone. 


To  study  of  the  performance  of  the  MEMS  device,  the  directional  cues  in  terms  of  the  mIPD  are 
characterized  with  respect  to  the  sound  frequency  and  the  sound  azimuth.  In  Figure  7,  it  can  be  seen  that 
mIPD  can  be  amplified  over  a  wide  frequency  range.  By  using  a  scanning  laser  vibrometer  (Polytec 
MSA400),  the  rocking  and  bending  modes  (see  the  mode  shapes  shown  in  Figure  7)  are  obtained  to  be  at 
10.3  kHz  and  29.5  kHz,  close  to  the  design  values.  Furthermore,  the  mIPD  is  seen  to  achieve  its 
maximum  value  and  experience  a  sign  change  at  the  rocking  natural  frequency,  while  at  the  bending 
natural  frequency,  the  mIPD  is  close  to  zero.  At  the  working  of  8  kHz,  the  experimental  value  of  mIPD  is 


clearly  a  linear  function  of  6  in  the  range  of  -3O°<0<3O°  (see  inset  of  Figure  7),  and  the  slope  of  the  linear 
range,  DS,  is  estimated  to  be  3.29  deg/deg,  which  is  18.3  times  of  the  initial  value  of  the  IPD 
(0.18deg/deg)  at  the  midline  without  amplification.  Such  a  performance  is  equivalent  to  that  can  be 
obtained  from  a  conventional  microphone  pair  with  a  separation  of  22.9  mm  (18  times  larger  than  that  of 
the  fly-ear  sized  sensor).  The  experimental  results  compare  well  with  the  simulation  results  obtained  on 
the  basis  of  a  2DOF  model. 

It  is  mentioned  that  the  damping  factors  of  the  device  are  determined  to  be  0.10  and  0.15,  which  are 
much  smaller  than  those  of  the  fly  ear.  Although  it  is  difficult  to  realize  the  damping  characteristics  of  the 
fly  ear  in  MEMS  devices,  smaller  damping  factors  render  a  higher  amplification  of  the  directional 
sensitivity  DS  (designed  value:  1.18  deg/deg  and  the  experimental  value:  3.29deg/deg),  which  can  help 
increase  the  signal-to-noise  ratio  for  sound  localization. 


Figure  7:  Directional  cue  mIPD  as  a  function  of  sound  frequency  for  the  fabricated  bio-inspired 
device.  The  mode  shapes  of  the  rocking  mode  (10.3  kHz)  and  the  bending  mode  (29.5  kHz)  obtained 
from  laser  Doppler  vibrometer  tests  are  also  shown.  The  simulation  results  are  obtained  by  using 
damping  factors  of  0.10  and  0.15.  The  initial  phase  difference  refers  to  the  acoustic  inputs.  At  the 
working  frequency  8  kHz,  the  measured  mIPD  has  a  linear  relationship  with  0  from  -30°  to  30°.  DS 
for  this  linear  azimuth  range  is  3.29  deg/deg,  18.3  times  the  initial  DS  at  the  midline  of  0.18  deg/deg. 


C.  Development  Miniature  Optical  Detection  and  Signal  Processing  System 

As  discussed  previously,  optical  detection  technique  has  been  used  to  interrogate  the  fly-ear 
inspired  sensors.  However,  even  with  a  miniature  size  of  the  sensor  itself,  the  large  scale  optical  system 
greatly  hinders  the  application  of  the  fly-ear  inspired  sensors.  In  this  work,  the  untapped  potential  of  a 
low  coherence  optical  sensor  system  is  explored  to  develop  a  Micro-Opto-Electro-Mechnical- System 
(MOEMS)  sensor  platform  that  is  capable  of  integrating  multiplexed  Fabry-Perot  (FP)  interferometer 
based  sensors.  A  schematic  of  the  proposed  differential  low  coherence  interferometer  based  sensing 
element  module  with  all  of  the  components  and  waveguide  integrated  on  a  single  MOEMS  chip  is  shown 
in  Figure  8.  Light  from  a  low  coherence  light  source  with  a  coherence  length  Lc  is  first  sent  to  the 
reference  interferometer  -  a  tunable  Fabry-Perot  interferometer,  and  then  via  a  waveguide  splitter  to  two 
combiners.  At  the  output  of  each  combiner,  the  light  beam  (Eini  or  Ein2)  is  sent  to  the  sensor  interferometer 
-  a  Fabry-Perot  interferometer  based  sensor.  A  Fabry-Perot  configuration  uses  two  mirrors  (a  partial 
mirror  and  a  second  mirror,  which  can  be  either  a  partial  or  a  complete  mirror)  to  form  a  cavity.  The 
reflected  light  waves  from  the  mirrors  (e.g.,  Eout2=  Eout2i+  Eout22)  will  have  an  optical  path  difference  (OPD) 
Ls  that  induces  a  phase  difference  <ps=koLs,  where  lu,  =2jt/A  is  the  free-space  wave  number  and  A  is  the 


Figure  8:  Schematic  of  low  coherence  interferometer  based  MOEMS  chip  integrated  with  biology- 
inspired  sensor. 

wavelength.  The  acoustic  pressure  induced  diaphragm  deflection  produces  a  phase  difference  change  A(j> 
in  cf)s.  The  reflected  light  from  each  Fabry-Perot  sensor  is  then  coupled  back  to  the  waveguide  and  sent  to 
a  photodetector.  The  on-chip  tunable  Fabry-Perot  interferometer,  which  has  an  initial  optical  path 
difference  Lr  (or  phase  difference  (pr=koLr),  acts  as  a  reference  interferometer.  When  it  is  path-matched  to 
the  Fabry-Perot  sensor  ( Lr  =LS)  and  the  coherence  length  LC«LT,  Ls ,  the  output  intensity  received  by  the 
detector  is 

lout  - 4  +  4  cos k0(L  -Lr)  =  4  +  4  cos (£  -</>r)  =  4  +  4  cos(A^) , 

where  A  and  B  are  the  constants  related  to  the  mirror  properties  of  the  FP  interferometer,  and  A(()  is  the 
differential  phase  change  between  the  sensing  interferometer  and  reference  interferometer.  Note  that  A(j)  is 
the  only  parameter  related  to  the  center  displacement  (X)  of  the  biology-inspired  diaphragm  and 
A(()=2koX.  Therefore,  the  pressure  sensitivity  (displacement/pressure)  of  the  diaphragm  can  be  amplified 
by  a  factor  of  2k0(107  times  at  /.=  1 300  nm).  Such  high  sensitivity  enables  the  detection  of  a  sound 
pressure  level  of  10  dB  (the  sound  level  of  calm  breathing).  In  addition,  owing  to  the  differentiation  of  the 
phase  signal  between  the  two  interferometers,  this  technique  has  immunity  to  wavelength  and  power 
fluctuation  induced  noise,  permits  a  short  effective  sensing  cavity  (several  pm),  and  yields  a  high 
resolution  (~  10'4  nm[5])  and  a  large  dynamic  range  (several  tens  of  wavelength).  This  is  the  first  attempt 
towards  developing  a  low  coherence  interferometer  based  MOEMS  detection  system. 

An  optical  Micro-Electro-Mechanical-System  (MEMS)  sensor  platform  was  developed  for  the 
miniature  fiber  optic  microphones.  The  key  components  of  this  platform  include  a  superluminescent 
diode  (SLED),  photodetector  chips,  and  a  MEMS  tunable  Fabry-Perot  (FP)  filter.  The  working  principle 
of  the  multi-functional  sensor  platform  is  based  on  a  differential  low  coherence  interferometer  consisting 
of  multiple  sensing  interferometers  (multiple  sensors)  and  a  single  reference  interferometer  (the  FP  filter). 
By  measuring  the  differential  phase  signal  between  the  reference  interferometer  and  each  of  the  reference 
interferometers,  this  platform  has  immunity  to  wavelength  and  power  fluctuation  induced  noise,  permits  a 
short  effective  sensing  cavity,  and  yields  a  high  resolution  and  a  large  dynamic  range. 

The  key  components  of  the  optical  MEMS  platform  are  shown  in  Figure  9.  The  SLED  can  emit  a 
broadband  light  with  a  center  wavelength  of  1310  nm  (coherent  length  35  pm)  and  a  power  of  120  pW  at 
100  mA.  The  tunable  FP  filter  is  formed  between  a  fiber  endface  and  a  curved  MEMS  mirror  driven  by  a 
comb  drive,  which  has  a  cavity  length  of  60  pm.  The  comb  drive  is  designed  to  have  a  resonant  frequency 
of  1.4  kHz.  The  photodiode  has  a  receiving  area  of  300  pm  in  diameter  and  a  responsivity  of  0.9  A/W  at 
1310  nm. 
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Figure  9:  Close-ups  of  key  components  of  the  optical  sensor  board  in  the  MEMS  platform:  (a) 
SLED;  (b)  Tunable  Fabry-Perot  filter;  (c)  photodiode. 

Finally,  we  have  integrated  the  MOEMS  on-chip  optical  system  with  a  wireless  transceiver 
system  to  realize  a  fully  functional,  inexpensive  and  portable  sensor  system.  The  final  prototype  of 
wireless  optical  signal  processing  system  is  shown  in  Figure  10. 


SLED  with  TEC  Photodiodes 


Figure  10:  Integrated  optical  detection  and  signal  processing  chip 


D.  Bio-Inspired  Sound  Localization  with  Fly-Ear  Inspired  Sensors 

A  bio-inspired  sound  localization  scheme  was  developed,  which  is  based  on  the  fly’s  unique 
localization/lateralization  scheme.  We  further  implemented  this  sound  localization  scheme  on  a  motorized 
control  platform  by  using  the  fly-ear  inspired  sensors  described  in  the  previous  section. 

Inspired  by  the  fly’s  unique  localization/lateralization  scheme  for  sound  source  localization,  a  bio¬ 
inspired  localization  scheme  is  developed  and  implemented  by  controlling  of  a  motorized  platform.  As 
shown  in  Figure  11,  this  platform  consists  of  two  motorized  one  dimensional  rotational  stages  (Newport 
URS75BPP),  with  one  stage  stacked  on  top  of  the  other.  The  bio-inspired  MEMS  acoustic  sensor 
discussed  in  the  previous  quarterly  report  was  mounted  on  the  tip  of  a  homemade  fixture.  A  speaker  (ESS 
AMT)  was  used  as  the  sound  source.  The  centers  of  both  stages  were  aligned  vertically  with  the  sensor  by 
manual  adjustment. 

The  bio-inspired  MEMS  directional  microphone  uses  the  low  coherence  fiber  optic  interferometer 
(LCFOI)  system  (see  Figure  6)  to  detect  the  displacement  of  the  two  membrane  centers.  The  two 
motorized  stages  are  connected  to  and  controlled  by  a  motion  controller  (Newport  ESP300).  The  outputs 
of  the  photo-detectors  are  connected  to  a  data  acquisition  board  (NI  USB-6259).  Both  the  motion 
controller  and  the  data  acquisition  board  are  connected  to  a  laptop  with  the  localization  algorithm 
implemented  in  a  Labview  program. 


Figure  11:  Control  platform  using  two  rotational  stages  and  integrated  with  a  bio-inspired 
MEMS  directional  microphone. 

First,  the  phase  difference  mIPD  of  the  MEMS  device  is  measured  at  various  azimuth  angles  when 
a  7kHz  pure  tone  (close  to  the  optimal  working  frequency  of  the  sensor)  is  played  through  the  speaker. 
The  mechanical  interaural  time  difference  (mlTD)  is  calculated  by  finding  the  maximum  of  the  cross¬ 
correlation  of  the  two  time  signals  for  the  two  membranes 

mlTD  =  arg  corr  ( Wt ,  Wr )  (11) 

which  is  converted  to  phase  difference  mIPD.  As  shown  in  Error!  Reference  source  not  found.,  mIPD  is  a 
linear  function  of  azimuth  in  the  range  of  -3O°<0<3O°.  By  approximating  their  relationship  in  the  linear 
range  by  a  straight  line,  the  directional  sensitivity  is  calibrated  as  1.32,  which  is  8.6  times  the  directional 
sensitivity  using  a  conventional  uncoupled  configuration.  It  should  be  noted  that  there  is  a  zero  offset 
(5.09°)  at  the  zero  azimuth  due  to  the  unsynchronized  data  acquisition  and  alignment  error.  In  short,  the 
phase  difference  in  the  linear  range  is  approximated  by  the  following  formula: 

mIPD  =  1. 32  x  (0-5.09°).  (12) 
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Figure  12:  Phase  difference  of  the  MEMS  sensor  as  a  function  of  azimuth  at  7kHz.  The  uncoupled 
case  corresponds  to  a  conventional  directional  microphone  setup  without  the  mechanical  coupling. 


Upon  receiving  two  signals  of  the  membranes,  the  control  program  first  calculates  the  phase 
difference  mlPD.  Then,  an  azimuth  estimate  is  obtained  by  using  the  following  formula: 


6  = 


mlPD 

1.32 


+  5.09°. 


(13) 


If  the  estimate  is  beyond  the  linear  range  (e.g.  45°),  the  control  program  moves  the  sensor  toward 
the  end  position  of  the  linear  range  (e.g.  30°).  Once  the  source  is  within  the  linear  range,  the  sensor  will 
be  moved  to  that  specific  location  so  that  it  will  point  to  the  source.  This  scheme  is  essentially  the  same 
as  the  fly’s  localization/lateralization  scheme. 

As  a  demonstration  of  the  bio-inspired  localization  scheme,  the  sensor  is  placed  at  80°  initially 
(speaker  stage  reading  of  0°  and  sensor  stage  reading  of  80°),  as  shown  in  Figure  13.  As  the  sensor 
pinpoints  the  source  (both  stage  readings  are  0°),  the  speaker  stage  starts  to  move  randomly  in  either 
direction.  However,  the  control  program  is  able  to  detect  the  random  motion  of  the  source  and  move  the 
sensor  stage  in  the  opposite  direction  so  that  the  sensor  is  maintained  to  be  able  to  still  pinpoint  the 
source. 
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Figure  13:  Real  time  speaker  and  sensor  position  readings  from  the  motorized  stages.  At  about  70 
interactions,  a  command  is  sent  to  the  speaker  stage  to  move  randomly  in  either  direction  and  the 
movement  step  increases  with  time. 

In  Figure  14,  the  actual  position  of  the  source  relative  to  the  sensor  is  compared  to  that  of  the 
estimate  value  from  the  localization  system.  In  the  beginning  stage,  the  source  is  out  of  the  linear  range, 
which  is  equivalent  to  30°  for  the  sensor.  Once  the  source  is  within  the  linear  range,  the  sensor  is  able  to 
locate  the  source  and  track  its  position  even  when  the  source  is  moving  in  a  random  fashion.  This 
experiment  successfully  demonstrates  the  application  of  fly-ear  localization/lateralization  scheme  in  one 
dimensional  sound  source  localization  and  tracking. 


Figure  14:  Time  history  of  actual  and  estimated  positions  of  the  sound  source. 


E.  Investigation  of  Optimal  Control  in  Fly’s  Localization  Scheme 

The  parasitic  fly  Ormia  ochracea  has  been  found  to  have  superacute  directional  hearing,  although 
its  distance  between  auditory  organs  is  only  half  millimeter.  In  addition,  the  fly  employs  a  unique 
localization/lateralization  scheme  to  locate  its  host  crickets:  the  fly  estimates  the  true  position  of  the 
sound  source  in  the  linear  range  of  azimuth  up  to  20-30deg,  whereas  it  only  makes  a  left/right  decision 
when  the  source  is  beyond  this  linear  range. 

We  previously  showed  that  this  scheme  may  relate  to  the  characteristics  of  its  directional 
sensitivity  (DS)  of  its  mechanically  coupled  eardrums.  At  its  working  frequency  5kHz,  which  is  the 
calling  song  frequency  of  male  crickets,  DS  is  not  only  flat  in  the  vicinity  of  midline  (azimuth  up  to 
30deg),  but  also  higher  than  what  is  achievable  at  any  other  frequencies.  Based  on  this,  an  argument  is 
made  that  the  fly  ear  represents  a  nature  designed  optimal  structure. 

The  goal  of  the  study  in  this  section  is  to  investigate  the  fly  ear  phenomena  another  angle,  from 
the  perspective  of  optimal  control.  By  modeling  in  ID  (azimuth),  the  optimal  control  problem  is  first 
described,  and  its  analytical  solution  is  obtained.  Then  efforts  are  made  to  explore  the  relationship 
between  the  fly’s  localization/lateralization  scheme  and  the  phase  portrait  of  the  optimal  control  solution. 

Description  of  optimal  control  problem 

•  Model  (bearing-only  ID  control  problem) 

Assuming  the  state  variables  are  the  bearing  angle  and  the  rotational  speed:  x1  =d,x2  =  6,  the 
state  equations  can  be  written  as: 

Xl  =  X, 

■  .  ,  (14) 

X2  =  U  -  ,UX2 

where  u  is  the  control  moment,  jlix2  is  the  dragging  force  that  is  assume  to  be  proportional  to  the 
rotational  speed. 

Another  assumption  is  on  the  control  moment,  which  is  subject  to  the  constraint: 

\u\  <  1. 


(15) 


This  constraint  limits  the  maximum  control  the  fly  can  apply.  As  a  result,  the  maximal  speed  is  equal 
to  1  /  fj. 

In  matrix  form,  (14)  can  be  written  as: 


X  =  AX  +  Bw, 


where 


X  =  J 
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Since  the  controllability  matrix  has  full  rank,  this  system  is  fully  controllable. 

•  Solution  of  state  equation 
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Taking  the  inverse  Laplace  transform  of  (20),  we  can  obtain  the  transition  matrix: 
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At 

e  = 
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Therefore,  the  state  at  any  time  t  is 

X(t)  =  eA'X(0)  +  j'  eA[,~T)Bu  (r]dr 
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If  the  control  u  is  constant  in  the  interval  ^t,t  +  AT ],  the  state  at  the  end  of  the  interval  is  related  to  the 
state  in  the  beginning  of  the  interval  by: 


where 


X(?  +  Ar)  =  AArX(?)  +  BAV 
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(23) 
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(25) 
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The  above  equation  is  helpful  for  the  numerical  implementation  in  Matlab. 

•  Performance  measure 

The  control  problem  here  is  to  move  from  any  initial  state  to  the  final  rest  state,  x  =  0  .  To  get  the 
optimal  control,  we  need  a  general  form  of  performance  cost: 

j(u )  -Jl'M-wi]  dt,  (26) 

which  is  a  weighted  sum  of  elapsed  time  and  control  effort  (or  consumed  fuel).  The  goal  is  to  find  the 
optimal  control  that  minimizes  the  performance  cost  J.  If  X  is  equal  to  zero,  J  reduces  to  control  effort 
only;  When  X  is  zero,  the  goal  is  simply  to  minimize  the  travel  time. 

For  the  optimal  control  problem  formulated  above,  an  analytical  solution  can  be  obtained  by 
using  Pontryagin’s  minimum  principle161,  as  follows. 


Solution  of  the  optimal  control  problem 

The  Hamiltonian  can  be  written  as: 

H  =  A  +  |w(*)|  +  pl(t)x2(t)  +  p2(t)[u(t)~  /ix2(t)  j 
=  |w(7)|  +  p2(t)u(t)  +  A  +  p{  {t)x2  (t)~  pp2  (t)x2  (7) 

where  px ,  p 2  are  the  co-state  variables. 

The  co- state  equations  are: 

■  *  /  x  dH  _ 

Pi  (0  =  -— =  0 

OXj 


The  solutions  of  the  co-state  are: 


Pi  {*)  =  -—  = -Pi  {*)  + Pi  {t)P  =  0 

cl\\ 


Pi  (0 = ci 


Pi(t)  =  ^  +  C2ef" 

For  the  optimal  control  u  [t ) ,  we  have  the  following  relationship  from  the  Pontryagin’s  minimum 
principle: 

|m*(?)|  +  p2*  (t)u*  (t)  <\u(t)\  +  p2*(t)u(t). 

We  only  need  to  choose  the  control  to  minimize  the  LHS  of  (30). 

Depending  on  the  value  of  p2,  the  control  u  can  be  take  one  of  the  three  values: 

[-1  P*i>\ 

u  (7)  = .  o  -l  <  p*2  <  l 

l  p2  <  -l 


Apparently,  the  trajectory  can  not  end  with  a  segment  withy*  =0.  Also  notice  that  p*  in  (29)  is  a 

♦ 

monotonically  changing  function  of  time  t,  which  means  u  can  not  change  abruptly  from  1  to  -1,  or  vice 
versa.  Therefore,  candidates  for  the  control  history  are: 

(1)  End  with  a  u  =  -1  segment:  j-lj,  jo  —  —  l} , {l  — >  0  — >  —  1  j- ; 

(2)  End  with  a  u  =  1  segment:  {l} , {0  —  l},{— 1  — *  0  — *  l}. 

However,  {-l}  and  jo— '-ljcan  be  considered  as  a  portion  or  a  special  case  ofjl  — 0  — -lj,  and 
jlj  and  |0  — 5 *  ljcan  be  considered  as  a  portion  of  j-1  — *  0  — *  l}.  Therefore,  we  only  need  to  consider 

two  control  scenarios:  {1-0—1}  and  {- 1  — >  0-1}  .  Starting  from  the  final  zero  state,  we  calculate 
the  state  set  at  which  the  control  makes  the  switch. 

•  Case  1:  j— 1  —  0  —  l} 

Assume  the  control  u  switches  from  -1  to  0  at  t\,  and  from  0  to  1  at  t2,  and  the  state  reaches  the  zero  state 
at  t{. 

(1)  Segment  A  <t  <tf:  u  =  1 


From  (22),  we  have 


xl(tf^  =  xl(t2)  +  —x2(t2)  1-e ^  — - 
X2  (tf  )  =  x2{t2)  +  -  [l  -  e^tf-h) '  = 


1  1 

—  +  ~ 
p  /d 


{tf  t2)  + 


1  =  Q 

P1 


Solving  (32)  gives: 


(^)  iStf  ^  p2 . 


(t2)  =  -\  1-e^ 
M  - 


Arll-e* 


Therefore,  the  state  lies  on  this  curve: 


*1  (*2  )  =  - "T  ln  [!  -  t*X2  (*2  )]  -  -  • *2  (*2  ) 


where  xl  {t2 )  >  0 and  x,  (t2 )  <  0 ,  since  p  >  0  and  tf  >  t2 . 
(2)  Segment  tl  <t  <t2:  u  =  0 


From  (22),  we  have: 

Xi  (t2  )  =  xl(t1)  +  -x2  (t, )  [l  -  e~4h~h)  ] 

X2  (^2  )  =  X2  (^1  )  ^ 

At  ti,  w*(^)<  0,  )  =  1 


(35) 


(36) 


^  +  c1e'H'  =1, 


H(fj )  =  A  +  Pi  (0*2  M~^x2  )  =  h  +  clx2  {tx)-px2  (/j )  =  0. 

At  ?2,  U  {t2  )  >  0  ,  /?2  (?2  )  = 

—  +  c7e/"2  =-l, 

H(?2 )  =  A  +  A  (f2 )x2  (?2 )  +  {*2 )  =A  +  cxx 2  (t2)  +  /.ix2  (t2 )  =  0. 

From  (36)  and  (38),  we  can  obtain  the  difference  t2-t{. 

1  ,  c,  +  u 

t2-tx  =  —  In  — - . 

it  c  j  -  u 

The  constant  Ci  can  be  solved  from  (37): 

A 

c\  ~  P  7~\  ■ 


x2(tx) 


1  .  ,  Z-  LI 

t2  -tx  =  —In  1 - x. 


Substitute  Ci  in  (40)  using  (41): 


Substituting  (42)  into  (35),  we  can  express  the  state  at  h  by  the  state  at  t\. 

L  (*  „  l+\  2x2  (0 


Xx(t2)  =  xx(tx)- 
X2{t2)  =  X2(tx)- 


A-2/jx2(tx) 


PV  2/  2  V '^-2^(0 

Apparently,  [xx  (t2 )  x2  [t2 )  j  is  on  the  trajectory  of  (34): 

/u  jU  A-2jUX2[tx) 

Where  xx  {tx )  >  Oand  x2  (tx )  <  0. 

•  Case  2:  jl  — »  0  — »  - 1  j 

Assume  the  control  u  switches  from  1  to  0  at  t\,  and  from  0  to  -1  at  t2,  and  the  state  reaches  the  zero  state 
at  t{. 

(1)  Segments,  <t  <tf:  u  =  - 1 
From  (22),  we  have 


Xi  (tf )  =  X,  (f, )  +  -  x,  (f2 )  [  1  -  e^'f "'2)  1  +  -4  -~{tf  - 12 )  = 


/A  LI 


x2  (tf )  =  x2  (t2  -  -  [l  -  '  = 


■  (45) 


The  switch  from  u  =  0  to  u  =  - 1  occurs  on  this  curve: 


(46) 


*1  {h )  =  -4ln  [1  +  //x2  (t2 )]  -jx2  (f2 ), 

where  xl  (t2 )  <  0 and  x2  (t2  )>(). 

(2)  Segment  tx<t  <t2:  u  =  0 


The  switch  from  u  =  1  to  u  =  0  occurs  on  this  curve: 

xM^^xMAln 

where  xx  (tt )  <  Oand  x2  (ij  )>(). 


1  + 


ukx2  ( tx ) 

A  +  2//X,  (f,) 


(47) 


Results  and  discussion 

Based  on  the  solutions  obtained  in  the  previous  section,  the  whole  state  can  be  split  into  four 
regions,  as  shown  in  Figure  15.  Starting  from  an  initial  state  point  A|  in  Region  I  (or  A2  in  Region  III), 
the  control  is  u  =  —1  (or  u  =  1 ),  till  the  state  reaches  the  first  switching  interface  OBi  (or  OB2),  after 
which  the  control  is  turned  off,  i.e.  u  =  0 ,  in  Region  II  (or  Region  IV).  When  the  state  approaches  the 
second  switching  interface  OCi  (or  OC2),  the  control  is  again  turned  on,  u  =  1  (or  u  =  -1 ).  Finally,  it 
reaches  the  final  zero  state  O  along  OCi  (or  OC2).  According  to  the  optimal  control  theory,  since 
Ai-*Bi->Ci-»0  is  the  optimal  trajectory  starting  from  initial  state  at  Ai,  Bl— »C]— »0  is  also  optimal  if 
the  trajectory  is  initiated  at  B[. 


Figure  15:  Exemplary  optimal  trajectories  (A  =  1,  //  =  0.01),  in  which  the  initial  state  is  at  A] 

(Xj  =  90, x2  =0)  or  A2  ( Xj  =  -90,  X,  =  0),  and  the  final  state  is  at  the  origin  O. 

Figure  16  shows  the  time  history  of  the  trajectory  Ai-»B]-»Ci-»0  in  Figure  15.  Starting 
fromXj  =  90, X,  =  0,  a  maximum  control  is  applied  (u  =  -l)  to  accelerate  toward  the  final  position 

X,  =  0.  After  coasting  in  the  middle  segment  of  the  trajectory,  the  system  is  controlled  by  a  deceleration 
force  u  =  1  to  reach  the  final  position  with  zero  velocity. 


(a) 


(a) 


Figure  16:  Time  history  of  the  trajectory 
A^Bj-fCi-fO  in  Figure  15  (A  =  l,ju  =  0.01). 
(a)  control,  (b)  azimuth,  and  (c)  velocity. 


Figure  17:  Time  history  of  the  trajectory 
starting  from  x]  =  90,  X1  =  0.  (a)  control,  (b) 

azimuth,  and  (c)  velocity.  The  parameters  of  the 
performance  measure  are  changed  to 

X  =  100,^  =  0,01. 


The  length  of  the  coasting  segment  varies  depending  on  the  relative  priority  of  time  and  control 
effort  in  the  performance  measure  in  (26),  which  is  controlled  by  the  parameter  ~k.  Increasing  ~k  from  1  to 
100  means  it  is  more  of  a  priority  to  reduce  the  travel  time.  The  time  history  in  Figure  17  shows  that  the 
total  travel  time  is  reduced  from  22  for  X=1  to  19  for  A.=100.  However,  this  reduction  is  at  the  expense  of 
increasing  control  effort.  Excluding  the  coasting  time,  the  time  for  u=±  1  instead  increases  from  1 1  for 
\=\  to  19  for  X=100. 

Revisit  of  localization/lateralization  scheme  of  the  fly  ear 

In  the  phonotactic  experiment  on  the  fly[4],  it  is  found  that  the  fly’s  turning  speed  is  a  sigmoid 
function  of  the  speaker  azimuth,  as  shown  in  Figure  18.  The  fly  employs  a  localization/lateralization 
scheme  to  accurately  locate  its  host  crickets,  i.e.  the  fly  estimates  the  true  position  of  the  source  when  the 
source  is  within  the  linear  range  (azimuth  up  to  20-30°),  and  otherwise  it  only  makes  a  left  or  right 
decision  and  turns  a  maximal  angle  toward  the  source. 

This  unique  scheme  is  related  to  the  fly’s  sensing  mechanism.  By  looking  at  the  directional 
sensitivity  (DS,  defined  as  the  first  derivative  of  mechanical  interaural  phase  difference  mIPD  with 
respect  to  azimuth  0),  as  shown  in  Error!  Reference  source  not  found.,  we  can  find  out  that  in  the  vicinity 
of  midline  (-3O°<0<3O°),  DS  is  not  only  very  flat  at  the  fly’s  working  frequency  5kHz,  but  it  is  higher 
than  what  is  achievable  at  any  other  frequency.  In  other  words,  if  we  use  a  straight  line  to  approximate 
mIPD  in  the  azimuth  range  -3O°<0<3O°,  the  linearity  is  best  at  5kHz,  and  the  slope  is  maximal.  The 
localization/lateralization  scheme  takes  full  advantage  of  the  linear  range  so  that  the  fly  can  accurately 
pinpoint  the  crickets. 


Speaker  azimuth  (degrees  relative  to  midline) 

Figure  18:  Fly’s  turning  size  measured  in  1.2  seconds  (equivalent  to  turn  speed)  as  function  of 
speaker  azimuth  (Mason  et  al,  2001). 

In  Figure  19,  the  phase  portrait  is  obtained  by  using  A  =  100,  //  =  0.1,  which  means  the  maximum  speed 
achievable  is  1  /  ^  =  10.  The  resulting  velocity-azimuth  relationship  is  very  similar  to  the  sigmoid 
relationship  in  Figure  18,  except  the  different  orientations  of  the  Z-shape.  This  is  simply  due  to  the  fact 
that  the  reference  is  the  fly  in  the  phonotactic  experiment,  while  the  reference  is  the  fixed  speaker  in  the 
current  simulation.  In  both  cases,  the  head  front  ends  up  facing  the  speaker. 


Figure  19:  Phase  portrait  for  the  ID  system  with  optimal  control  for  initial  state  at  Ai(90,-10)  or 
A2(-90,10).  The  parameters  for  the  optimal  measures  are  A  =  100,  ju  =  0.1. 


Inspired  by  the  fly’s  localization/lateralization  scheme,  we  try  to  investigate  how  to  design  a  control  that 
is  close  to  the  optimal  control.  Assuming  the  full  state  is  observable,  one  example  of  closed-loop  control 
is  as  follows: 

r-1  6  >30 


u  —  i 


-0.5x  — -0.5,mx,  -30  <6  <30 

30  2 


1  6  <  -30 

where  the  control  is  still  subject  to  the  constraint  \u\  <  1 . 


(48) 


Figure  20  compares  the  results  of  the  above  defined  closed-loop  control  with  the  optimal  control.  For  the 
closed-loop  control,  the  state  overshoots  the  target:  when  xi  reaches  0  for  the  first  time  at  t=9.7,  x2  is 
equal  to  -6.38.  Although  it  settles  downs  in  the  end,  it  takes  68.7  time  to  get  to  the  zero  state,  much  longer 

than  13  by  the  optimal  control.  Also,  the  control  effort  (defined  as  ^\u\dt)  is  29.6%  more  (14.58  as 

opposed  to  1 1.25). 

Note  that  in  the  beginning,  the  feedback  is  based  on  the  position  xi  only  so  that  only  the  azimuth  needs  to 
be  estimated.  In  this  case,  there  is  no  additional  damping  introduced  by  the  feedback;  the  real  part  of  the 
closed-loop  system  is  just  -u.  Therefore,  the  system  will  take  longer  to  settle  down. 


(a)  (b) 


(c)  (d) 


Figure  20:  Comparison  of  the  optimal  control  (blue  solid  line)  and  the  closed-loop  control  (red 
dashed  line)  in  (35). 
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Abstract  (Andrew  Lisiewski) 

A  micro-scale  sound  localization  sensor  is  developed  and  studied  in  this  thesis  to  address  the  fundamental 
challenge  of  miniaturizing  sound  localization  systems.  When  miniaturizing  a  microphone  array  there  is  a 
critical  size  limitation  at  which  the  array  will  be  unable  to  localize  the  sound  source  in  a  discernible 
manner.  However,  a  solution  to  this  dilemma  came  about  when  studying  the  hearing  mechanisms  of  a 
particular  fly,  known  as  Ormia  ochracea.  Background  research  into  the  hearing  mechanisms  of  the  fly 
found  that  it  can  accurately  locate  a  sound  source  even  though  its  eardrums  are  separated  by  a  distance  of 
only  500  pm.  The  fly’s  exceptional  directional  hearing  capability  has  been  linked  to  a  distinct  mechanical 
coupling  between  its  two  eardrums,  which  helps  amplify  minute  directional  cues.  Inspired  by  the 
remarkable  hearing  capabilities  of  the  fly’s  micro-scale  ear,  researchers  have  sought  to  develop  micro¬ 
scale  sensors  to  mimic  the  fly’s  ear.  One  limitation  of  simply  imitating  the  fly’s  ear  is  that  the  fly  is  only 
capable  of  localizing  a  sound  source  in  one  dimension.  In  this  thesis  work,  the  knowledge  gained  from 
understanding  the  fly  ear  mechanism  is  applied  to  achieve  the  goal  of  developing  a  micro-scale  sound 


localization  sensor  capable  of  sound  source  localization  in  two  dimensions.  In  this  thesis,  for  the  first 
time,  micro-scale  fly-ear  inspired  sensor  devices  employing  three  or  four  coupled  membranes  have  been 
designed.  Reduced-order  models  have  been  developed  to  achieve  a  fundamental  understanding  of  the 
performance  of  each  sensor  design.  Furthermore,  a  micro-scale  sensor  device  incorporating  three 
mechanically  coupled  membranes  arranged  in  an  equilateral  triangular  configuration  has  been 
successfully  developed.  Experimental  study  of  the  sensor  device  incorporated  with  a  low  coherence  fiber 
optic  interferometric  detection  system  has  suggested  that  the  micro-scale  fly-ear  inspired  sensor  can 
achieve  a  much  improved  performance  in  terms  of  phase  differences  and  directional  cues  when  compared 
to  a  similar  sized  microphone  array  constructed  with  separate  microphones.  In  addition,  localization 
techniques  have  been  developed  to  best  use  the  fly-ear  inspired  sound  localization  sensors.  Future  work 
is  suggested  to  incorporate  this  sensor  system  with  a  fully  autonomous  robot  to  improve  robot  homing 
and  navigation. 

Andrew  Lisiewski  has  received  his  Master’s  degree  in  Spring  2011. 

Abstract(Haiiun  Liu) 

Microphone  arrays  have  been  widely  used  in  sound  source  localization  for  many  applications  including 
hearing  aid  devices,  robot  navigation,  and  underwater  sensor  networks.  In  order  to  locate  the  sound  in  a 
discernible  manner,  the  separation  between  microphones  needs  to  be  greater  than  a  critical  distance, 
which  poses  a  fundamental  challenge  for  the  miniaturization  of  directional  microphones.  In  nature,  the 
auditory  organs  of  the  parasitoid  fly  Ormia  ochracea  are  forcibly  set  close  to  each  other  (about  500  pm), 
but  it  can  localize  its  cricket  hosts  with  a  resolution  of  as  small  as  2°.  The  key  to  this  remarkable 
directional  hearing  is  that  the  two  tympana  (ear  drums)  are  coupled  by  a  cuticular  bridge  that  is  pivoted 
about  its  center.  Due  to  such  mechanical  coupling  mechanism,  the  time  difference  is  amplified  by  more 
than  30  times,  and  the  best  intensity  difference  can  reach  10  dB.  This  innovative  solution  in  nature  can 
inspire  one  to  find  alternative  approaches  to  developing  miniature  directional  microphones. 

The  overall  goal  of  this  dissertation  work  is  to  unravel  the  underlying  physics  of  the  fly  ear  hearing 
mechanisms,  and  to  apply  this  understanding  to  develop  novel  bio-inspired  miniature  directional 
microphones.  In  the  first  step  of  this  proposed  work,  the  fly  ear’s  localization  performance  will  be 
investigated  based  on  an  equivalent  two-degrees-of-freedom  model.  Secondly,  a  novel  bio-inspired 
directional  microphone  design  with  mechanically  coupled  diaphragms  will  be  proposed  to  capture  the 
essential  dynamics  of  the  fly  ear.  For  this  design,  a  non-dimensionalized  continuum  mechanics  model  will 
be  developed.  Parametric  studies  will  be  carried  out  to  explore  how  the  key  normalized  parameters  affect 
the  performance  of  this  directional  microphone.  Such  non-dimensionalization  makes  the  analysis 
independent  of  the  physical  size,  large  or  small.  Next,  this  mechanics  model  will  be  used  to  guide  the 
development  of  large-scale  proof-of-concept  microphones,  which  will  be  experimentally  studied  by  using 
a  low-coherence  fiber  optic  interferometric  detection  system.  Finally,  miniaturization  efforts  will  be  made 
to  develop  microscale  devices  by  using  micro-fabrication  techniques  along  with  an  on-chip  integrated 
optical  detection  system. 

It  is  expected  that  this  dissertation  work  will  provide  physical  insights  into  the  fly  ear  mechanism  and 
sound  localization  physics.  The  developed  bio-inspired  miniature  directional  microphone  would 
significantly  improve  the  localization  performance  while  bringing  down  the  size  limit  as  imposed  on  the 
traditional  sensors.  As  such,  it  would  pave  the  way  to  new  applications  that  have  constrained  space  but 
require  high  sensitivity. 

Haijun  Liu  is  expected  to  defend  his  Ph.D.  dissertation  in  late  Spring  2012. 


Yuxiang  Liu  has  worked  on  the  optical  detection  part  of  the  project,  whose  dissertation  was  on  a  different 
topic.  He  received  his  Ph.D.  in  Fall  2010. 
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